Anticipated Synchronization in a Biologically Plausible Model of Neuronal Motifs 
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Two identical autonomous dynamical systems coupled in a master-slave configuration can exhibit 
anticipated synchronization (AS) if the slave also receives a delayed negative self- feedback. Recently, 
AS was shown to occur in systems of simplified neuron models, requiring the coupling of the neuronal 
membrane potential with its delayed value. However, this coupling has no obvious biological corre- 
late. Here we propose a canonical neuronal microcircuit with standard chemical synapses, where the 
delayed inhibition is provided by an interneuron. In this biologically plausible scenario, a smooth 
transition from delayed synchronization (DS) to AS typically occurs when the inhibitory synaptic 
conductance is increased. The phenomenon is shown to be robust when model parameters are varied 
within physiological range. Since the DS-AS transition amounts to an inversion in the timing of the 
pre- and post-synaptic spikes, our results could have a bearing on spike-timing-dependent-plasticity 
models. 

PACS numbers: 87.18.Sn, 87.19.11, 87.19.1m 



I. INTRODUCTION 

Synchronization of nonlinear systems has been exten- 
sively studied on a large variety of physical and biological 
systems. Synchronization of oscillators goes back to the 
work by Huygens, and in the past decades an increased 
interest in the topic of synchronization of chaotic systems 
appeared [lj. 

About a decade ago, Voss Q discovered a new scheme 
of synchronization that he called "anticipated synchro- 
nization" . He found that two identical dynamical sys- 
tems coupled in a master-slave configuration can exhibit 
this anticipated synchronization if the slave is subjected 
to a delayed self-feedback. One of the prototypical exam- 
ples proposed by Voss is described by the equations 
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f(x) is a function which defines the autonomous dynam- 
ical system. The solution y(t) = x(t + td), which char- 
acterizes the anticipated synchronization (AS), has been 
shown to be stable in a variety of scenarios, including 
theoretical studies of autonomous chaotic systems , 
inertial ratchets [f| , and delayed-coupled maps Q , as well 
as experimental observations in lasers @, H[ or electronic 
circuits [|| . 

More recently, AS was also shown to occur in a non- 
autonomous dynamical system, with FitzHugh-Nagumo 
models driven by white noise fl(il - [T2j ]. In these works, 
even when the model neurons were tuned to the excitable 
regime, the slave neuron was able to anticipate the spikes 
of the master neuron, working as a predictor . Though 
potentially interesting for neuroscience, it is not trivial to 
compare these theoretical results with real neuronal data. 
The main difficulty lies in requiring that the membrane 
potentials of the involved neurons be diffusively coupled. 
While a master-slave coupling of the membrane poten- 
tials could in principle be conceived by means of electri- 



cal synapses (via gap junctions) (l3j or ephaptic interac- 
tions [14j . no biophysical mechanism has been proposed 
to account for the delayed inhibitory self-coupling of the 
slave membrane potential. 

In the brain, the vast majority of neurons are cou- 
pled via chemical synapses, which can be excitatory or 
inhibitory. In both cases, the coupling is directional and 
highly nonlinear, typically requiring a suprathrcshold ac- 
tivation (e.g. a spike) of the pre-synaptic neuron to 
trigger the release of neurotransmitters. These neuro- 
transmitters then need to diffuse through the synaptic 
cleft and bind to receptors in the membrane of the post- 
synaptic neuron. Binding leads to the opening of spe- 
cific channels, allowing ionic currents to change the post- 
synaptic membrane potential 13(. This means that not 



only the membrane potentials are not directly coupled, 
but the synapses themselves are dynamical systems. 

Here we propose to bridge this gap, investigating 
whether AS can occur in biophysically plausible model 
neurons coupled via chemical synapses. The model is 
described in section |TTJ In section IIHI we describe our 
results, showing that AS can indeed occur in "physio- 
logical regions" of parameter space. Finally, section IIVI 
brings our concluding remarks and briefly discusses the 
significance of our findings for neuroscience, as well as 
perspectives of future work. 



II. MODEL 

A. Neuronal motifs 

We start by mimicking the original master-slave cir- 
cuit of eqs. (JTJ) with a unidirectional excitatory chemical 
synapse (M — > S in Fig.QJi). In a scenario with standard 
biophysical models, the inhibitory feedback we propose 
is given by an interneuron (I) driven by the slave neu- 
ron, which projects back an inhibitory chemical synapse 
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to the slave neuron (see Fig. [TJl) . So the time-delayed 
negative feedback is accounted for by chemical inhibition 
which impinges on the slave neuron some time after it 
has spiked, simply because synapses have characteristic 
time scales. Such inhibitory feedback loop is one of the 
most canonical neuronal microcircuits found to play sev- 
eral im por tant roles, for instance, in the spinal cord [l5j . 
cortex 15 1, thalamus [l6l Il7j and nuclei involved with 
song production in the bird brain [l8[ . For simplicity, we 
will henceforth refer to the 3-neuron motif of Fig. [TJi as 
a Master-Slave-Interneuron (MSI) system. 



(a) 



(b) 




FIG. 1. (a) (Color online) Three neurons coupled by chemical 
synapses in the master-slave-interneuron (MSI) configuration 
: excitatory AMPA synapses (with maximal conductance #a) 
couple master (M) to slave (S) and slave to interneuron (I), 
whereas an inhibitory GABAa synapse (with maximal con- 
ductance go) couples interneuron to slave, (b) Same as (a), 
except that all three neurons of the MSI circuit receive exci- 
tatory (NMDA) synapses from a driver neuron (D). 

As we show in section IIIII below, whether or not the 
MSI circuit can exhibit AS depends, among other fac- 
tors, on the excitability of the three neurons. In the 
MSI, this is controlled by a constant applied current (see 
section UlI A[) . To test the robustness of the results (and 
at the same time improve the realism and complexity of 
the model), in section IlIIB I we study the four-neuron mo- 
tif depicted in Fig. [TJd, where the excitability of the MSI 
network is chemically modulated via synapses project- 
ing from a global driver (D). From now on, we refer to 
the 4-neuron motif as a Driver-Mastcr-Slave- Interneuron 
(DMSI) microcircuit. 



B. Model neurons 

In the above networks, each node is described by a 
Hodgkin-Huxley (HH) model neuron [l9j], consisting of 
four coupled ordinary differential equations associated to 



the membrane potential V and the ionic currents flowing 
across the axonal membrane corresponding to the Na, K 
and leakage currents. The gating variables for sodium 
are h and m and for the potassium is n. The equations 
read H: 



C m-jr = G Na m 3 h{E 
at 



dx 
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Na-V) + G K n\E K -V) 
V) + I + J2lsyn (2) 
-P*{V)x, (3) 



where x £ {h, m, n}, C m = 9ir fjF is the membrane ca- 
pacitance of a 30 x 30 x it /im 2 equipotential patch of 
membrane [2~oT | , / is a constant current which sets the neu- 
ron excitability and ^ I syn accounts for the interaction 
with other neurons. The reversal potentials are E^ a = 
115 mV, E K = -12 mV and V rest ^10.6 mV, which 
correspond to maximal conductances Gnu — 1080-7T mS, 
Gk = 3247T mS and G m = 2.7ir mS, respectively. The 
voltage dependent rate constants in the Hodgkin-Huxley 
model have the form: 
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Note that all voltages are expressed relative to the resting 
potential of the model at / = [20]. 

According to Rinzel and Miller [2l|, in the absence 
of synaptic currents the only attractor of the system of 
equations (H])-© for / < 177.13 pA is a stable fixed point, 
which loses stability via a subcritical Hopf bifurcation at 
I ~ 276.51 pA. For 177.13 pA < I < 276.51 pA, the 
stable fix point coexists with a stable limit cycle. 



C. Synaptic coupling 

AMPA (A) and GABAa (G) arc the fast excitatory 
and inhibitory synapses in our model [see Fig. [TJl] . Fol- 
lowing Destexhe et al [22|, the fraction (i — A,G) 
of bound (i.e. open) synaptic receptors is modelled by a 
first-order kinetic dynamics: 



dt 



_ r (0)_ 



(0 



(10) 



where and (3i are rate constants and [T] is the neuro- 
transmitter concentration in the synaptic cleft. For sim- 
plicity, we assume [T] to be an instantaneous function of 
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the prc-synaptic potential V pre : 

^ v ^= l+e Jv::-v v)/KpV (id 

where T max = 1 mM -1 is the maximal value of [T], K p = 
5 mV gives the steepness of the sigmoid and V p = 62 mV 
sets the value at which the function is half-activated [22| . 
The synaptic current at a given synapse is given by 

jW =g i r^(V-E i ), (12) 

where V is the postsynaptic potential, gt the maximal 
conductance and Ei the reversal potential. We use Ea — 
60 mV and E G = -20 mV. 

The values of the rate constants a a, Pa, eta, and Pg 
are known to depend on a number of different factors 
and vary significantly [23U2H ] . To exemplify some of our 
results, we initially fix some parameters, which are set 
to the values of Table 1 unless otherwise stated (sec- 
tion flTTXT]). Then we allow these parameters (as well 
as the synaptic conductances) to vary within the phys- 
iological range when exploring different synchronization 
regimes (see section lmA~2l and HlTB|) . 

The slow excitatory synapse is NMDA (N) and its 
synaptic current is given by: 

=g N B(Vy N \V -E N ), (13) 

where En = 60 mV. The dynamics of the variable r^ N ' 
is similar to eq. (fTU)) with ajv = 0.072 mM _1 ms _1 and 
p N = 0.0066 his -1 . The magnesium block of the NMDA 
receptor channel can be modeled as a function of postsy- 
naptic voltage V: 

B (y^ = 2 7 e (-0.062V)[Mg 2 +] o /3.57 ' 

where [Mg 2+ ] Q = 1 mM is the physiological extracellular 
magnesium concentration. 

In what follows, we will drop the neurotransmitter su- 
perscripts A, G and TV from the synaptic variables r and 
I. Instead we use double subscripts to denote the re- 
ferred pre- and postsynaptic neurons. For instance, the 
synaptic current in the slave neuron due to the interneu- 
ron (the only inhibitory synapse in our models) will be 
denoted as lis, and so forth. 





MSI 


DMSI 


otA (mM i ms x ) 


1.1 


1.1 


Pa (ms" 1 ) 


0.19 


0.19 


olg (mM _1 ms _1 ) 


5.0 


5.0 


Pa (ms" 1 ) 


0.30 


0.60 


a N (rnM^W) 




0.072 


Pn (ms" 1 ) 




0.0066 


9A (nS) 


10 


10 


7(pA) 


280 


160 



TABLE I. Standard values employed in the model. See text 
for details. 



III. RESULTS 
A. Master-Slave-Interneuron circuits 

1. Three dynamical regimes 

Initially, we describe results for the scenario where all 
neurons receive a constant current / ^ 280 pA. This 
corresponds to a situation in which the fixed points are 
unstable and, when isolated, they spike periodically. All 
other parameters are as in Table 1. For different sets of 
inhibitory conductance values gc our system can exhibit 
three different behaviors. To characterize them, we define 
tf 1 as the time the membrane potential of the master 
neuron is at its maximal value in the i-th cycle (i.e. its 
i-th spike time), and tf as the spike time of the slave 
neuron which is nearest to tf 1 . 

The delay t; is defined as the difference (see Fig. [5]): 

n = if - tf. (15) 

Initial conditions were randomly chosen for each com- 
puted time series. When Tj converges to a constant value 
r, a phase- locked regime is reached [27|. If r < ("mas- 
ter neuron spikes first" ) we say that the system exhibits 
delayed synchronization (DS) [Fig.[2Ja)]. If r > ("slave 
neuron spikes first"), we say that anticipated synchro- 
nization (AS) occurs [Fig. [2fb)]. If r does not converge 
to a fixed value, the system is in a phase drift (PD) 
regime ■ The extent to which the AS regime can be le- 
gitimately considered "anticipated" in a periodic system 
will be discussed below. 

In Figure[3]we show examples of time series in the three 
different regimes (DS, AS and PD). The different panels 
correspond to the membrane potential, fraction of acti- 
vated receptors for each synapse, and synaptic current in 
the slave neuron. For a relatively small value of the in- 
hibitory coupling [go = 20 nS, Fig.^a)] the slave neuron 
lags behind the master, characterizing DS. In Fig. [3£b), 
we observe that by increasing the value of the inhibitory 
coupling (go = 40 nS) we reach an AS regime. Finally, 
for strong enough inhibition [gc — 60 nS, Fig. |3{c)] the 
PD regime ensues. 

In the DS and AS regimes the master and slave neu- 
rons spike at the same frequency. However, when the 
system reaches the PD regime the mean firing rate of the 
slave neuron becomes higher than that of the master. 
The counterintuitive result shown in Fig. |4ja) emerges: 
the mean firing rate of the slave neuron increases while 
increasing the conductance of the inhibitory synapse pro- 
jected from the interneuron. For the particular combina- 
tion of parameters used in Fig. IHa), the transition turns 
out to be reentrant, i.e., the system returns to the DS 
regime for sufficiently strong inhibition (a more detailed 
exploration of parameter space will be presented below) . 
Figure H^b) shows the return map of the inter-spike in- 
terval of the slave, which forms a closed curve (touching 
the trivial single-point return map of the master). This 
is consistent with a quasi-periodic phase-drift regime. 
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FIG. 2. (Color online) Membrane potential V as a function 
of time for an external current / = 280 pA in the master (M), 
slave (S), and interneuron (I) neurons. The plot illustrates 
two regimes: (a) <?g = 20 nS leads to delayed synchronization 
(DS), where r < 0, and (b) go = 40 nS leads to anticipated 
synchronization (AS), where r > 0. Other parameters as in 
Table 1. 



Note that in this simple scenario go plays an analo- 
gous role to that of K in Eq. [T] for which AS is stable 
only when K > K c (eventually with reentrances) [28j |. 
Moreover, the behavior of the synaptic current in the 
slave neuron is particularly revealing: in the DS regime 
[Fig. Efa)], it has a positive peak prior to the slave spike, 
which drives the firing in the slave neuron. In the AS 
regime [Fig. E^b)], however, there is no significant re- 
sulting current, except when the slave neuron is already 
suprathreshold. In this case, the current has essentially 
no effect upon the slave dynamics. This situation is sim- 
ilar to the stable anticipated solution of Eq. [TJ when the 
coupling term vanishes. 
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FIG. 3. (Color online) Time series of the membrane poten- 
tials (V), bound receptors (r) and synaptic currents (I), with 
model parameters as in Table 1. Note that the system is peri- 
odic in the DS and AS regimes [(a) and (b) respectively], but 
not in the PD regime (c). 



2. Scanning parameter space 

The dependence of the time delay r on go is shown in 
Fig. [5] for different values of the external current / and 
maximal excitatory conductance gA- Several features in 
those curves are worth emphasizing. First, unlike previ- 
ous studies on AS, where the anticipation time was hard- 
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FIG. 4. (Color online) (a) The mean firing rate of the slave 
(Fs) coincides with the mean firing rate of the master (Fm) 
for DS and AS regimes, but it is larger for PD. (b) In PD, the 
return map of the inter-spike interval of the slave is consistent 
with a quasi-periodic system (the pink star shows the return 
map of the master). 

wired via the delay parameter td [see eq.©], in our case 
the anticipation time r is a result of the dynamics. Note 
that go (the parameter varied in Fig. [5]) does not change 
the time scales of the synaptic dynamical variables (r), 
only the synaptic strength. 

Secondly, r varies smoothly with go- This continu- 
ity somehow allows us to interpret r > as a legiti- 
mately anticipated regime. The reasoning is as follows. 
For go = 0, we simply have a master-slave configuration 
in which the two neurons spike periodically. Due to the 
excitatory coupling, the slave's spike is always closer to 
the master's spike which preceded it than to the mas- 
ter's spike which succeeded it [as in e.g. Fig. ©(a)]. 
Moreover, the time difference is approximately 1.5 ms, 
which is comparable to the characteristic times of the 
synapse. In that case, despite the formal ambiguity im- 
plicit in the periodicity of the time series, the dynamical 
regime is usually understood as "delayed synchroniza- 
tion" . We interpret it in the following sense: the system 
is phase-locked at a phase difference with a well defined 
sign [27}. Increasing go, the time difference between the 
master's and the slave's spikes eventually changes sign 
[as in e.g. Fig. ©(b)]. Even though the ambiguity in 
principle remains, there is no reason why we should not 
call this regime "anticipated synchronization" (again a 
phase-locked regime, but with a phase difference of op- 
posite sign). In fact, we have not found any parameter 
change which would take the model from the situation 
in Fig. ©(a) to that of Fig. ©(b) by gradually increas- 
ing the lag of the slave spike until it approached the next 
master spike. If that ever happened, r would change 
discontinuously (by its definition). Therefore, the term 
"anticipated synchronization" by no means implies vio- 
lation of causality and should just be interpreted with 
caution. As we will discuss in section IIV1 the relative 



timing between pre- and postsynaptic neurons turns out 
to be extremely relevant for real neurons. 

Third, it is interesting to note that the largest antic- 
ipation time can be longer (up to 3 ms, i.e. about 20% 
of the interspike interval) than the largest time for the 
delayed synchronization (« 1.5 ms). If one increases go 
further in an attempt to obtain even larger values of r, 
however, the system undergoes a bifurcation to a regime 
with phase drift [which marks the end of the curves in 
Fig. EJ. 
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FIG. 5. (Color online) Dependence of the time delay r with 
the maximal conductance go for different values of the applied 
current / and gA ■ The end of each curve (stars) marks the 
critical value of go, above which the system changes from AS 
to PD. 

The number of parameters in our model is very large. 
The number of dynamical regimes which a system of 
coupled nonlinear oscillators can present is also very 
large, most notably p/g-subharmonic locking structured 
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FIG. 6. (Color online) Delay r (right bar) in the (fiu, <?g) pro- 
jection of parameter space: DS (blue, right), AS (red, middle) 
and PD (white, left — meaning that no stationary value of r 
was found). 



P A = 0.19 ms" 1 




10 20 30 40 50 

S A (nS) 



10 20 30 40 50 



10 20 30 40 50 

Sa(i>S) 



FIG. 7. (Color online) Delay r (right bar) in the (<7a,<?g) 
projection of parameter space for different combinations of 
Pa and pa- From left to right we have respectively PD, AS 
and DS regimes, as in Fig. [6j 
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FIG. 8. (Color online) Delay r (right bar) in the (<M,Sg) 
projection of parameter space for different values of I. PD, 
AS and DS regimes as in Fig. [5] 



in Arnold tongues [29( . These occur in our model, but not 
in the parameter region we are considering. In this con- 
text, an attempt to map all the dynamical possibilities in 
parameter space would be extremely difficult and, most 
important, unproductive for our purposes. We there- 
fore focus on addressing the main question of this paper, 
which is whether or not AS can be stable in a biophysi- 
cally plausible model. 

In Fig. [5] we display a two-dimensional projection of 
the phase diagram of our model. We employ the val- 
ues in Table U except for gA , which is varied along the 
horizontal axis. Note that each black curve with circles 
in Fig. [5] corresponds to a different vertical cut of Fig. [51 
along which go changes. We observe that the three differ- 
ent regimes arc distributed in large continuous regions, 
having a clear transition between them. Moreover the 
transition from the DS to the AS phase can be well ap- 
proximated by a linear relation golSA ~ 3.5 in a large 
portion of the diagram. 

Linearity, however, breaks down as parameters are fur- 
ther varied. This can be seen e.g. in Fig. [7j which dis- 
plays the same projection as Fig. [51 but for different com- 
binations of Pg and Pa- We observe that AS remains 
stable in a finite region of the parameter space, and this 
region increases as excitatory synapses become faster. 

Figure [5J suggests that larger values of the input cur- 
rent I eventually lead to a transition from AS from DS. 
This effect is better depicted in Fig. [51 where the DS re- 
gion increases in size as / (and therefore the firing rate) 
increases. Figurcs[5|b)-(d) also show that the system can 
exhibit reentrant transitions as go is varied. Most impor- 
tantly, however, is that Figs. [7J and [5] show that there is 
always an AS region in parameter space, as synaptic and 
intrinsic parameters are varied. 

As we will discuss in section ITVl the possibility of con- 
trolling the transition between AS and DS is in prin- 
ciple extremely appealing to the study of plasticity in 



neuroscience. However, in a biological network, the in- 
put current would not be exactly constant, but rather be 
modulated by other neurons. In the following, we test the 
robustness of AS in this more involved scenario, therefore 
moving one step ahead in biological plausibility. 



B. Driver-Master-Slave-Interneuron circuits 

Let us consider the MSI circuit under a constant input 
current I = 160 pA. This is below the Hopf bifurca- 
tion [ll[ , i.e. none of the three neurons spikes tonically. 
Their activity will now be controlled by the driver neu- 
ron (D), which projects excitatory synapses onto the MSI 
circuit [see Fig. [IJb)]. We chose to replace the constant 
input current by a slowly varying current, so that the 
synapses projecting from the driver neuron are of the 
NMD A type (see section [II]). The driver neuron receives 
a current Id = 280 pA, so it spikes tonically. All remain- 
ing parameters are as in the second column of Table 1. 
The interest in this case is to verify whether AS holds 
when the excitability of the MSI circuit is modulated by 
a non-stationary current. 

As shown in Fig. [51 we found in this new scenario a 
similar route from DS to AS, and then the PD regime 
(compare with Fig. [Jj. Note that the characteristic time 
(/3jv = 6.6s~ 1 )for the unbinding of the NMD A receptors 
is about ten times larger than the inter-spike interval of 
the driver neuron (which spikes at » 67 Hz). As a conse- 
quence, r^M, fus, tdi are kept at nearly constant values 
(with variations of « 10% around a mean value — data 
not shown). The variations in the NMDA synaptic cur- 
rent are also small, which in principle should make the 
system behave in an apparently similar way to the pre- 
vious MSI circuit. However, these small variations are 
important enough to increase the AS domain in parame- 
ter space, in some cases even eliminating the PD region 
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FIG. 9. (Color online) DMSI circuit (see Fig. [Da). Delay r 
(right bar) in the (<7a,<?g) projection of parameter space for 
different combinations of /3a and /3g- PD, AS and DS regimes 
as in Fig. [6] 
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FIG. 10. (Color online) DMSI circuit (see Fig.[Tp). (a) The 
mean firing rate of the slave (Fs) coincides with the mean 
firing rate of the master (Jm) for DS and AS regimes, but it 
is larger for PD. (b) In PD, the return map of the inter-spike 
interval of the slave is consistent with a quasi-periodic system. 



(see e.g. Fig. |9]for fla = 0.30 ms" 1 ). Therefore, at least 
in this case, the use of more biological plausible parame- 
ters docs not destroy AS, but rather enhances it. 

In fact, the three regions in the MSI diagrams seem 
to retain their main features in the DMSI circuit. When 
PD occurs, for example, the slave again spikes faster than 
the master [Fig. 110( a)]. like in the MSI circuit [compare 
with Fig. Ufa)]. Another signature of the robustness of 
the PD phase against the replacement of a constant by 
a slowly-varying synaptic current appears in the return 
map shown in Fig. [TDT b) . It can be seen that it has the 
same structure of its three-neuron counterpart shown in 



Fig. lib). 



IV. CONCLUDING REMARKS 

In summary, we have shown that a biologically plausi- 
ble model of a 3-neuron (MSI) motif can exhibit an at- 
tractor in phase space where anticipated synchronization 
is stable. The transition from the DS to the AS regime is 
a smooth function of the synaptic conductances. Typi- 
cally, a further increase in the inhibitory conductance go 
leads to a second transition from AS to PD, a quasipcri- 
odic regime in which the slave firing frequency is larger 
than that of the master. 

We have varied synaptic decay rates (/3), synaptic con- 
ductances (g) as well as input currents (/) within well ac- 
cepted physiological ranges [23M2fj|. In all the scenarios 
there is always a continuous region in parameter space 
where AS is stable. Replacing the constant current by 
a global periodic driver (arguably a more realistic sit- 
uation), we obtain a model of a 4-neuron (DMSI) mo- 
tif which exhibits the same three regions of the simpler 
model. The synaptic rise constants (a) were also var- 
ied, but have a lesser effect on the transitions among the 
different regimes (data not shown). Therefore the phe- 
nomenon seems to be robust at the microcircuit scale. 

It is important to emphasize that our AS results differ 
from those obtained from eq. ([l} at a fundamental level. 
In our model the delayed feedback that leads to AS is 
given by biologically plausible elements (an interneuron 
and chemical synapses). Hence, the anticipation time is 
not hard-wired in the dynamical equations, but rather 
emerges from the circuit dynamics. Moreover, the par- 
ticular circuit we study is a neuronal motif ubiquitously 
found in the brain pjJ-Ojl • We are unaware of other AS 
models in which every parameter has a clear biological 
interpretation. 

We believe that our results can be extremely rele- 
vant for modeling studies of synaptic plasticity. Recent 
decades have witnessed a growing literature on spike- 
timing dependent plasticity (STDP), which accounts 
for the enhancement or diminution of synaptic weight 
[long term potentiation (LTP) and long term depression 
(LTD), respectively] depending on the relative timing be- 
tween the spikes of the pre- and post-synaptic neurons 
(see e.g. j30l - [33 |). Experimental data strongly suggest 
that if the pre-synaptic neuron fires before (after) the 
post-synaptic neuron, the synapse between them will be 
strenghtened (weakened) [33, [34|. STDP is supposed 
to take place in a window of time differences between 
post- and pre-synaptic spikes in the order of ten millisec- 
onds, within which the delay and anticipation times of 
our models fall. Since the DS-AS transition amounts to 
an inversion in the timing of the pre- and post-synaptic 
spikes, then by appropriately controlling this effect one 
could dynamically toggle between synaptic strengthening 
and weakening. This could be potentially linked with 
modeling of large-scale ascending feedback modulation 
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from reward systems. 

Our results, therefore, offer a number of possibilities 
for further investigation. Including effects from microcir- 
cuit dynamics (such as the ones we have presented here) 
in models of synaptic plasticity is a natural next step, one 
which we are currently pursuing. Once we have verified 
AS in a biologically plausible model, one could consider 
using simplified models [35|,[36[ (e.g. by replacing the HH 
equations and/or the synaptic kinetics) and the influence 
of noise 0,0. We are also investigating whether the 
structure of the phase diagram can be qualitatively repro- 
duced via a phase-response-curve analysis [37], HH of the 
neuronal motifs studied here. Results will be published 



elsewhere. 
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